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^ \ Previous studies have demonstrated that putatively single nitrogen-type Wolf- 

Rayet stars (WN stars) without known companions are X-ray sources. However, 
almost all WN star X-ray detections so far have been of earlier WN2 - WN6 

^^ \ spectral subtypes. Later WN7 - WN9 subtypes (also known as WNL stars) have 

proved more difficult to detect, an important exception being WR 79a (WN9ha). 
We present here new X-ray detections of the WNL stars WR 16 (WN8h) and 
WR 78 (WN7h). These new results, when combined with previous detections, 
demonstrate that X-ray emission is present in WN stars across the full range of 
spectral types, including later WNL stars. The two WN8 stars observed to date 
(WR 16 and WR 40) show unusually low X-ray luminosities (L^) compared to 
other WN stars, and it is noteworthy that they also have the lowest terminal 
wind speeds (vqo). Existing X-ray detections of about a dozen WN stars reveal 
a trend of increasing L^. with wind luminosity \jwind = |Mv^, suggesting that 
wind kinetic energy may play a key role in establishing X-ray luminosity levels 
in WN stars. 
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Subject headings: stars: individual (WR 16; WR 78; WR 79a) — stars: Wolf- 
Rayet — X-rays: stars 



1. Introduction 

The detection of X-ray emission from Wolf-Rayet (WR) stars was one of the many no- 
table discoveries of the Einstein Observatory (Seward et al. 1979; Pollock 1987). Subsequent 
studies have broadened our knowledge of the X-ray properties of WR stars but the under- 
lying physical mechanisms ultimately responsible for their X-ray emission are still not well 
understood. Most current theoretical models attribute the X-ray emission of WR stars (and 
0-type stars) to thermal plasma produced in shocks that form in their powerful supersonic 
winds. WR stars have typical terminal wind speeds Vq^ ~ 1000 - 5000 km s~^ and mass-loss 
rates M ~ 10~^ Mq yr^^. It has also been suggested that WR stars might produce non- 
thermal X-rays via such processes as inverse-Compton scattering, synchrotron emission, or 
nonthermal bremsstrahlung (Pollock 1987; Chen & White 1991). However, most WR stars 
detected so far show line-dominated X-ray spectra that can be acceptably fitted by thermal 
models. But, nonthermal emission has not yet been ruled out in a few cases (e.g. WR 142; 
Sokal et al. 2010). 

Several different shock-related processes could in principle lead to thermal X-ray emis- 
sion in WR stars. First, soft X-rays (kT < 1 keV) may be produced in shocks embedded 
in the outflowing wind that form as a result of the line-driven instability (Lucy & White 
1980; Lucy 1982; Feldmeier et al. 1997; Gayley & Owocki 1995). This mechanism may be 
responsible for the soft X-ray emission detected in some 0-type stars (Cassinelli et al. 2001; 
Kahn et al. 2001). But, all WR stars detected so far show a harder component (kT > 2 
keV) in addition to a soft component that is usually present. The hard component is not ex- 
pected if the emission arises solely from line-driven instability shocks. Second, a sufficiently 
strong magnetic field could confine the wind and channel it into two opposing streams which 
collide near the magnetic equator, forming a magnetically-confined wind shock (MCWS). 
This mechanism is capable of producing hard X-rays (kT ^ 2 keV) and may be relevant 
to X-ray production in young magnetic OB stars such as 9^ Ori C (Babel & Montmerle 
1997b; Gagne et al. 2005) and a Ori E (Skinner et al. 2008). However, very strong surface 
magnetic fields would be needed to confine the more powerful winds of WR stars (ud-Doula 
& Owocki 2002; Skinner et al. 2010, hereafter SIO), and observational evidence for strong 
magnetic fields on WR stars is currently lacking. Thus, the question of whether the MCWS 
mechanism operates in WR stars remains open. Finally, in the case of WR + OB binary 
systems, hard X-rays (kT > 2 keV) can originate in the region between the two stars where 
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the interaction of their winds results in a colliding wind (CW) shock (Prilutskii & Usov 1976; 
Luo et al. 1990; Stevens et al. 1992; Usov 1992). Hot X-ray plasma that is thought to arise 
at least partially in a CW shock has now been detected in grating spectra of several X-ray 
bright WR + OB binaries (7^ Vel: Skinner et al. 2001; Schild et al. 2004; WR 140: Pollock 
et al. 2005; WR 147: Zhekov & Park 2010). 

X-ray spectra of WR stars are needed to distinguish between the above possibilities. 
Although good-quality X-ray grating spectra of a few X-ray bright WR+OB binaries have 
been obtained (as noted above), the study of putatively single WR stars as a class has been 
hampered by a general lack of X-ray spectral data. The observational situation has improved 
considerably during the past few years as a result of the completion of an exploratory X-ray 
survey of single WR stars with XMM-Newton and Chandra (SIO). Moderate-resolution CCD 
X-ray spectra are now available for several WN stars but attempts to detect carbon-rich WC 
stars have so far only yielded upper limits on their X-ray luminosities with typical constraints 
of log Lx < 31.0 ergs s^^ (Skinner et al. 2006; see also Fig. 9 of SIO). 

A surprising result of the initial survey was that all detected WN stars showed evidence 
of hot plasma {kThot ~ 2 keV) in their X-ray spectra. This result is at odds with X-ray 
production via the line-driven instability mechanism, for which only soft emission is expected. 
It is thus apparent that other mechanisms are needed to fully explain the origin of X-ray 
emission in WN stars and, within the current theoretical framework, either the MCWS or 
CW shock processes could be responsible. However, there is no compelling justification to 
invoke the CW interpretation because the WN stars in the initial survey sample were selected 
on the basis of lack of evidence for binarity. The presumption that these stars are indeed 
single is of course subject to the usual caveat that close companions can escape detection. 

In the initial survey, XMM-Newton and Chandra data for the following detected WN 
stars were analyzed: WR 2 (WN2), WR 18 (WN4), WR 20b (WN6h), WR 24 (WN6ha), 
WR 134 (WN6), WR 136 (WN6h), and WR 79a (WN9ha). In addition, archival ROSAT 
data were analyzed for WR 16 (WN8h) and WR 78 (WN7h). Neither of these two WNL 
stars G showed significant X-ray emission (see also Sanders et al. 1985), although a marginal 
{2a) excess was noted in hard-band ROSAT PSPCcI images of WR 78. We present here 



^Nitrogen-rich WR stars with spectral types in the range WN7 - WN9 are referred to as late WN stars, 
or WNL stars. Those with spectral types of WN2 - WN5 are called early WN stars or WNE stars. WN6 
stars are intermediate and may be classified as either WNE or WNL (Crowther 2007). A "h" suffix in the 
spectral type denotes the presence of hydrogen emission lines in the spectrum. 

^The i? 05^ T Position Sensitive Proportional Counter (PSPC) provided energy coverage in the wO.l - 
2.35 keV energy range but lacked sensitivity at higher energies E > 2.35 keV. PSPC hard-band images cover 
the range E w 0.5 - 2 keV. In contrast, XMM-Newton provides coverage to higher energies and we use the 
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more sensitive XMM-Newton observations of WR 16 and WR 78 that were obtained after 
the resuhs of our initial WN survey were pubhshed (SIO). These observations are important 
because they provide additional information on the X-ray properties of WNL subtypes, 
which comprise ~24% of all WR stars listed in the Vllth Catalog of Wolf-Rayet Stars (van 
der Hucht 2001). WNL stars are believed to be the least chemically-evolved WR stars and 
some 0-type stars probably enter the WR phase as WNL stars (Crowther & Bohannan 
1997). WNL stars typically have higher mass-loss rates but lower terminal wind speeds than 
WNE stars (Table 2 of Crowther 2007). We report here a faint X-ray detection of WR 16 
(WN8h) and bright emission including a hot plasma component from WR 78 (WN7h). 



2. Observations 

The XMM-Newton observations are listed in Table 1. The general properties of the 
WNL stars discussed here are summarized in Table 2 and their X-ray properties are given in 
Table 3. The targets were selected from the Vllth WR Catalog (van der Hucht 2001) based 
on proximity, low visual extinction, and lack of evidence for binarity. In addition to the 
new detections of WR 16 and WR 78, we also include data from our previous XMM-Newton 
observation (SIO) of the WNL star WR 79a for comparison. 

Data were acquired with the European Photon Imaging Camera (EPIC) in full-window 
mode. EPIC provides charge-coupled device (CCD) imaging spectroscopy from the pn cam- 
era (Striider et al. 2001) and two nearly identical MOS cameras (MOSl and M0S2; Turner 
et al. 2001). The EPIC cameras provide energy coverage in the range E ^ 0.2 - 15 keV 
with energy resolution E/AE ^ 20 - 50. The MOS cameras provide the best on-axis angular 
resolution with FWHM ^ 4.3" at 1.5 keV. 

Data were reduced using the XMM-Newton Science Analysis System (SAS vers. 11.0) 
using standard procedures including the filtering of raw event data to select good event 
patterns and removal of data within time intervals of high background radiation. Spectra 
and light curves were extracted from a circular region of radius r = 15" (~68% encircled 
energy) centered on the source, thus capturing most of the source counts while reducing 
background counts. Background analysis was conducted on circular source-free regions near 
the source. The SAS tasks rmfgen and arfgen were used to generate source-specific RMFs 
and ARFs for spectral analysis. The data were analyzed using the HEASOFT Xanadu 
software package. 



energy range E = 2 - 8 keV to refer to the XMM-Newton hard-band. 
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3. Results 

3.1. WR 16 (= V396 Car = HD 86161) 

This WN8h star has the lowest terminal wind speed (vqo ~ 630 - 740 km s~^) of any 
WN star in our exploratory X-ray survey. The star is surrounded by a cocoon of molecular 
gas which was probably ejected during a previous evolutionary phase (Marston et al. 1999). 
Optical variability is present but its cause is not fully understood (Antokhin et al. 1995). An 
early study attributed the variability to an unseen low-mass companion (Moffat & Niemela 
1982) but this interpretation was later called into question (Manfroid et al. 1987). The 
stellar distance is somewhat uncertain with published values ranging from 2.37 - 4.36 kpc 
(Crowther et al. 1995; Nugis Sz Lamers 2000; van der Hucht 2001, Hamann et al. 2006) and 
we adopt here an intermediate value d ~ 3.6 kpc. 

Visual inspection of energy-filtered EPIC images showed no obvious signs of an X-ray 
source at the stellar optical position. However, more in-depth image analysis revealed a 
weak emission excess above the local background level in both the pn and MOS images. We 
extracted counts (0.3 - 8 keV) within a circular region of radius 15" (68% encircled energy) 
centered on the star in each EPIC detector, and background counts in different source-free 
regions near the star. After background subtraction we obtained 31 ±6 net pn counts in 
30.0 ks of usable exposure and 24±5 net MOSl+2 (sum of MOSl and M0S2) in 66.0 ks of 
usable exposure (33.0 ks per MOS). The faint emission is clearly visible in the contoured 
MOSl-l-2 image in Figure 1. Although the statistical significance of the excess is low {p^2a), 
its presence in all three EPIC detectors at an offset of <1" from the HST GSC optical position 
(Table 3) strengthens the conclusion that weak X-ray emission from WR 16 is detected. 

The emission is too faint for spectral and timing analysis. However, we used the PIMMS 
U simulator to estimate the X-ray luminosity of WR 16 based on the count rates. The pn 
rate of 1.03 c ks~^ gives an unabsorbed luminosity log Lx(0.3 - 8 keV) = 31.14 [30.77 - 31.30] 
ergs s~^ for an assumed distance d = 3.6 kpc, where the range in square brackets reflects 
the range of published distances d = 2.37 - 4.36 kpc. A PIMMS estimate based on the MOS 
count rate gives a nearly identical unabsorbed value log Lx(0.3 - 8 keV) = 31.17 [30.80 - 
31.33] ergs s~^. For the underlying spectrum, we assumed an absorbed two-temperature 
optically thin thermal plasma (apec) with kTi = 0.6 keV and kT2 = 3.5 keV, typical of 
detected WN stars (SIO). We adopted an absorption column density Nh = 3.7 x 10^^ cm~^ 
based on Ay = 1.67 (Table 1) and the Gorenstein (1975) conversion from Ay to Nh. The 



^For information on PIMMS (Portable Interactive Multi-Mission Simulator) see 



littp://cxc. harvard.edu/ciao/alielp/pimms. html . 
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slightly different conversion of Vuong et al. (2003) decreases the upper limit on Lx by 0.05 
dex. 

The above Lx is the lowest of any WN star detected so far, but the derived value 
is dependent on the assumed spectral model and absorption. Most WN stars show X-ray 
absorption in excess of that expected from Ay (SIO; see also the results for WR 78 in Sec. 3.2 
below). If that is the case for WR 16, then its unabsorbed Lx could be larger than the value 
inferred above. In this regard, we note that the pn counts for WR 16 are distributed roughly 
equally between the soft (0.3 - 2 keV) and hard (2-8 keV) bands. But, XSPEC simulations 
show that if the absorption were as low as expected based on Ay, only about ~15% - 20% 
of the counts would emerge in the hard band. This suggests that excess absorption above 
that expected from Ay is present in WR 16, but a higher signal-to-noise spectrum would be 
needed to confirm this. For those WN stars where sufficient counts are available for spectral 
analysis, the best-fit Nh values are typically ?»2.0 - 2.5 times larger than expected from 
estimates based on Ay (Table 4 of SIO). If we assume that is the case for WR 16 then its 
unabsorbed Lx determined from PIMMS is RiO.lO - 0.15 dex larger than the values quoted 
above. 



3.2. WR 78 (= V919 Sco = HD 151932) 

This WN7h star has been included in numerous optical studies and there is fairly good 
agreement in published stellar parameters when different studies are compared (Table 2). 
Optical variability has been reported (Vreux et al. 1987; Bratschi & Blecha 1996) but no 
periodicity that would hint of a companion has been found. 

XMM-Newton provides the first unambiguous detection of WR 78 (Fig. 2). The source 
is quite hard as is apparent from its median photon energy E50 = 2.47 keV (Table 3). Clear 
evidence for high-temperature plasma is provided by the detection of the Fe Ka emission 
line complex at 6.67 keV (Fig. 3). This line forms at a characteristic temperature T ~ 40 
MK. Using pn events in the 0.3 - 8 keV range, the Kolmogorov-Smirnov (KS) test gives a 
probability of constant count rate Peons = 0.12 and a x^ test applied to the background- 
subtracted pn light curve binned at 1000 s intervals gives Peons = 0.97. The KS test operates 
on all detected events within the source extraction region, and thus includes both background 
and source events. The lower value of Peons based on the KS test may be due to background 
effects. But, neither the KS test nor x^ analysis suggests statistically significant variability 
for WR 78 during the 26.26 ks (7.3 hours) of usable exposure. Likewise, the other WN stars 
in the survey have failed to show any clear X-ray variability in observations spanning f5i3 - 
10 hours (SIO). 
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Acceptable fits of tlie CCD X-ray spectra of WN stars detected in tlie initial survey were 
obtained using a two-temperature (2T) optically thin plasma model consisting of cool (kTi) 
and hot (kT2) plasma components plus a single absorption component (Nh). This model 
is expressed symbolically as NH*(kTi + kT2), where parentheses denote that both plasma 
components are viewed through the same absorption column. Interestingly, this model gives 
an unacceptable fit of the WR 78 X-ray spectrum as does a simpler IT model of the form 
NH*kTi. These single-absorption models fit the spectrum above 2 keV (including the Fe 
Ka line) reasonably well, but generally underestimate the observed flux at energies below 
1 keV. A similarly poor flt below 1 keV occurs for other single-absorption models such as 
a plane-parallel shock model {vpshock in XSPEC). In order to obtain an acceptable flt for 
WR 78, different absorption components were required for the cool and hot components, i.e. 
a model of the form NH,i*kTi + NH,2*kT2. 

The acceptable two-absorber model is summarized in Table 4. The inferred plasma 
temperatures kTi = 0.64 [0.48 - 0.88; 90% conf.] keV and kT2 = 2.27 [1.63 - 2.90] keV are 
not substantially different from the WN9h star WR 79a (Sec. 3.3) when 90% confldence 
ranges are taken into account, but the value of kT2 is at the low end of the range for all WN 
stars observed so far (SIO). Most of the X-ray flux and emission measure (as gauged by the 
value of norm in Table 4) of WR 78 is associated with the hotter plasma component, which 
is evidently viewed through much higher absorption than the cool component. The separate 
contributions of the cool and hot components are shown in the bottom panel of Figure 3. 
The best-flt absorption values of Nh,i = 1.38 [1.20 - 1.60] x 10^^ cm^^ and Nh,2 = 6.55 [4.10 
- 10.7] X 10^^ cm~^ are substantially larger than expected on the basis of visual extinction 
esimates, since Ay = 1.55 mag (Table 2) corresponds to Nh = 3.4 x 10^^ cm~^ using the 
conversion of Gorenstein (1975). The two-absorber model for WR 78 gives an unabsorbed 
X-ray luminosity log La.(0.3 - 8 keV) = 32.84 ergs s~^ (d = 1.99 kpc; Table 4), which is 
at the high end of the range for detected WN stars (Fig. 4; SIO). However, this value is 
quite uncertain because of the large dereddening correction. All spectral models that we 
considered, including simplistic IT models, give values log L2;(0.3 - 8 keV) > 32.1 ergs s~^ 
which we consider to be a reliable lower limit. 

The most likely cause of the anomalous X-ray absorption is the WR wind. In that case, 
the higher absorption associated with the hot plasma component suggests that it originates 
deeper in the wind (closer to the star) where the wind density is higher. A similar conclusion 
was reached for the WN6h star WR 136 (SIO) and for 0-type stars studied by Waldron & 
Cassinelli (2007) and Naze (2009). In the ideal case of a smooth unclumped wind, the radius 
of optical depth unity (as measured from the star) decreases toward higher X-ray energies 
so harder photons are able to escape from smaller radii near the star. Softer X-ray photons 
produced far out in the wind suffer less wind absorption and can escape whereas any soft 



photons produced in deeper denser layers of the wind are absorbed and remain undetected. 
In the more complex case of an inhomogenous (clumped) wind, the above picture is not 
altered substantially if the clumps are optically thin to X-rays but could well be different if 
dense optically thick clumps are present. 



3.3. WR 79a (= HD 152408) 

This WN9ha star is the latest WN subtype detected in X-rays so far. It is discussed in 
more detail in SIO but for comparison we summarize its X-ray properties in Table 3 and give 
its best-fit spectral parameters in Figure 3. The EPIC pn spectrum of WR 79a is softer than 
that of WR 78 as evidenced by its lower median photon energy (Table 3) and it is apparent 
from Figure 3 that the X-ray absorption of WR 79a is less than that of WR 78. The pn 
spectra of both stars show the high-temperature S XV (2.46 keV) emission line but the very 
high temperature Fe Ka complex (6.67 keV) detected in WR 78 is not seen in WR 79a. 



4. Discussion 

4.1. What Stellar Parameters Control WN Star X-ray Luminosity? 

The fundamental stellar parameters that govern intrinsic X-ray luminosity levels in WR 
stars have not yet been identified. Good quality X-ray data or useful upper limits exist for 
only about a dozen putatively single WN stars, so robust statistics and stringent tests for 
correlations of L^ against stellar parameters are still limited by small sample size. This is 
in contrast to massive 0-type stars, from which WR stars evolve. X-ray data now exist for 
hundreds of 0-type stars and the empirical relation L^ ~ 10~'^Li)oi was already becoming 
apparent from a large sample of O and early B stars detected in the i? 05"^ T All- Sky Survey 
(Berghofer et al. 1997). 

Based on existing XMM-Newton and Chandra X-ray data, there is no convincing evi- 
dence for a similar L^ oc Lboi relation in putatively single WN stars. Likewise, no evidence 
for such a relation was found using ROSAT data (Wessolowski 1996). As we have already 
shown (SIO), WN stars with similar values of L^oi can have L^ values differing by an order 
of magnitude or more. Most detected WN stars, including WR 78 (Table 4), have L^ > 
10^'^Lhoi with typical values being L^ ~ 10^^'^Lf,o;, as shown in Figure 4. However there are a 
few WNL stars whose Lx/Lboi ratios are at least an order of magnitude lower. Most notable 
among these is WR 16 (WN8h), for which we obtain log Lx/Lboi = —8.2 ± 0.2 when taking 
the range of published distances and stellar luminosities into account (Table 2) . A similarly 
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low 'Lx/'Lboi ratio is also assured for WR 40 (WN8h) by virtue of its XMM non-detection 
with an X-ray upper limit log La,(0.5 - 10 keV) < 31.6 ergs s"^ at an assumed distance d = 
3 kpc (Gosset et al. 2005). 

Wind shock models of X-ray emission from massive stars predict that L^^ will depend 
on mass-loss parameters. In the line-driven instability shock model, \jx oc M" where a = 1 
for the radiative shock case, a = 2 for the adiabatic case, and a = 0.6 if the observed L^ 
oc Lhoi scaling relation in 0-type stars is to be reproduced (Owocki et al. 2011). From an 
observational standpoint, there is no clear dependence of L^, on M for WN stars. In particular 
WR 2 (WN2) has one of the lowest mass-loss rates tabulated for WN stars (log M = —5.4 
- —5.3 Mq yr~^; Nugis & Lamers 2000; Hamann et al. 2006) but its L^ is comparable to 
that of other WN stars having M values an order of magnitude larger. In addition, the two 
WN8 stars WR 16 and WR 40 are underluminous in X-rays by about an order of magnitude 
compared to other WN stars with similar mass-loss rates. Thus, mass-loss rate alone is not 
a reliable predictor of L^ in WN stars. 

There is one clue in the existing data that suggests L^ in WN stars may depend on 
terminal wind speed Vqo or some function thereof. Specifically, we note that the two observed 
WN8 stars WR 16 (voo ~ 630 - 740 km s^^) and WR 40 (voo ~ 650 - 910 km s^^) have the 
lowest terminal wind speeds of all WN stars (Hamann et al. 2006; see also Gosset al. 2005 
for WR 40). These two stars also have the lowest L^ values, even though only an upper 
limit is currently available for WR 40. The WNSh star WR 156 also has a similar low wind 
speed (voo ~ 660 km s~^; Hamann et al. 2006) but it has not yet been observed in X-rays. 
In the absence of a binary companion we expect that any X-ray detection of WR 156 would 
be quite faint. 

The apparent falloff in L^. with decreasing Vqo suggests that L^ may be related to wind 
momentum (= Mvoo) or wind luminosity L^md (= |Mv^). Using data in the initial WN X- 
ray survey (SIO), we already noted a trend of increasing L^ with L^ind- The new detections of 
WR 16 and WR 78 presented here conform to and strengthen this trend (Fig. 5). This trend 
was not seen in ROSAT WN-star data analyzed by Wesselowski (1996), perhaps because 
ROSAT had little sensitivity above ~2 keV where a significant fraction of the X-ray flux in 
WN stars emerges. Using the data in Figure 5, the generalized Kendall's tau test gives a 
correlation probability Pcorr = 0.93. A linear regression fit gives log L^ = 0.89(±0.47)log 
^wind + 26.41 (±0.55) where L^ is in units of ergs s~^ and L^j^^ in units of lO^^M© yr~^ 
X (km s~^)^. Thus, the regression fit gives a nearly linear dependence L^, oc L^^^j'^'^^ but 
the slope is quite uncertain due to scatter in the data. Factors that could contribute to this 
scatter include uncertainties in distances. X-ray luminosity dereddening corrections, WR 
star mass-loss rates (which are sensitive to clumping effects), and terminal wind speeds. 
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A notable outlier is the WN6I1 star WR 136, which appears to be underluminous in 
X-rays compared to other WN stars with comparable Lyjin^ values (Fig. 5) and to other WN 
stars with comparable stellar luminosities Lgtar (Fig. 4). Its Lgtar value is also low compared 
to other WNh stars (Fig. 6 of Hamann et al. 2006). This raises a question as to whether 
current distance estimates may be too low. The estimates range from d = 1.26 kpc (Hamann 
et al. 2006) to d = 1.82 kpc (Nugis & Lamers 2000), whereas our derived X-ray luminosity 
log L^ = 31.51 ergs s~^ is based on the Hipparcos distance d = 1.64 kpc (SIO). We note that 
there is also a large uncertainty in the L^ value of WR 136 because of poorly-constrained 
X-ray absorption (Nh) in spectral fits. Some statistically-acceptable spectral models require 
higher absorption and yield dereddened values as high as log L^ ~ 32.2 ergs s~^ (d = 1.64 
kpc). Such higher values, if correct, would shift the position of WR 136 upward in Figures 
4 and 5 (as shown by the L^, error bar), in better agreement with other WN stars. 

If the above L^ oc Luimd relation refiects a true physical dependence, then the obvious 
conclusion would be that a relatively constant fraction of the kinetic energy in WN star winds 
is being converted into X-ray energy. Interestingly, the MCWS model of Babel & Montmerle 
(1997a) predicts such a dependence. Specifically, they obtain L^ ^ j^^lhock^ ^^^ t^^y i^^fer 
Vshock ~ Voo in order to reproduce the X-ray luminosity of the AOp star IQ Aur from their 
model. Thus, there is some theoretical support for a L^. oc L^^jind relation in early-type stars 
with magnetic fields. Even so, some observational scatter in this relation would be expected 
due to differences in B-field strengths between stars (eq. [10] of Babel & Montmerle 1997a) 
and bearing in mind that collision of the two hemispherical wind components could in some 
cases occur at subterminal speeds. Although the existence of strong surface magnetic fields 
on single WR stars is not yet established, some WR + OB binaries do exhibit nonthermal 
radio emission which may be associated with weak magnetic fields above the surface (Abbott 
et al. 1986; Skinner et al. 1999; Dougherty & Williams 2000). Since the MCWS mechanism 
operates in single stars there is no need to postulate unseen companions in this interpretation. 

If close unseen companions are present in our sample of putatively single WN stars then 
hot plasma could be created via the colliding wind mechanism. If the WR wind momentum 
dominates that of the companion, then the CW shock would form at or near the companion 
surface. This could also give rise to a L^, oc L^md dependence, but considerable scatter 
would be expected since in this case the predicted L^ also depends on binary separation and 
the radius of the companion star (Usov 1992; Luo et al. 1990). Since we have no specific 
information on the existence of such companions, their stellar or wind properties, or their 
orbital parameters, a meaningful comparison of the observed L^ values with CW model 
predictions cannot yet be made. 
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4.2. High X-ray Plasma Temperatures in WN Stars 

All WN stars for which good-quality CCD X-ray spectra exist show a high-temperature 
plasma component (kT hot ~ 2 keV) and almost always a cooler component (kT coo/ < 1 keV). 
The cool component could be due to radiative instability wind shocks, but the origin of the 
hotter component is not yet known. Analysis of XMM EPIC and Chandra ACIS spectra of 
WN stars reveals a broad range of hot-component temperatures from kThot ~ 2 keV (WR 
1, WR 78) up to kThot f« 4 - 8 keV (WR 20b, WR 110, WR 134; SIO). 

Such high-temperature plasma could be produced by supersonic WR winds shocking 
onto denser downstream material, or the wind from the opposite hemisphere (as in the 
MCWS model), or the wind or surface of as yet undetected companions (CW shocks). These 
models all predict that the maximum shock temperature should scale proportionally with 
the pre-shock wind-speed as kT^™"^"^ ^ 1.96yuv^^, keV, where /j, is the mean atomic weight 
per particle and vi;^. is the wind speed in units of 1000 km s~^ (Babel & Montmerle 1997a; 
Luo et al. 1990;, Stevens et al. 1992). For solar abundances /j, ~ 0.6 but for metal-enriched 
WN star winds n ~ 1.06 (van der Hucht et al. 1986). For the WN star case of interest here 
one thus obtains kT^™"^^ ^ 2voo,ifc^ keV, assuming that the wind has reached terminal speed 
at the shock interface. 

Terminal wind speeds of X-ray detected WN stars lie in the range Voo ~ 700 - 2500 km 
s~^ (Table 2; SIO). Thus, the above shock models predict maximum shock temperatures 
'^'^ shock ~ 1 - 12 keV. The values of kThot determined from X-ray spectra lie comfortably 
within this range. However, there is no clear trend for a proportionality relation kThot oc v^, 
as one might naively expect from shock-model predictions. The WN star with the lowest 
Voo for which a good-quality CCD spectrum is available is the WN9ha star WR 79a (vqo ~ 
955 km s^"*^. Table 2) and its best-fit temperature is kThot = 2.66 [1.66 - 5.60; 90% conf.] 
keV (SIO). By comparison, the WN4 star WR 1 has a higher terminal speed (vqo ~ 1600 - 
2135 km s^^; Nugis & Lamers 2000; Ignace et al. 2001; Hamann et al. 2006) but a lower 
best-fit temperature kThot = 1-72 [1.36 - 2.21; 90% conf.] keV based on our analysis of XMM 
archive data (obsid 0552220101; see also Ignace et al. 2003a for analysis of an earlier XMM 
observation with a shorter exposure). 

Given the above results, it is not yet clear that terminal wind speed is the primary 
parameter governing hot-plasma component temperatures in WN stars. Other factors such 
as magnetic fields (if present) could play a role in plasma heating. Even if the hot-component 
plasma temperature is governed by Vqo, establishing such a dependence from the existing 
data will not be straightforward. As noted above, shock models predict that the maximum 
shock temperature kT^™"^^^ oc v^, but shocks will in general produce plasma over a range of 
temperatures. X-ray spectral fits of CCD spectra using simplistic 2T models only provide 
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an estimate of the average hot-component temperature kThot, which will in general be less 
than kT)j™^"^ . Furthermore, the inferred values of Whot from fits of moderate resolution 
CCD spectra can be uncertain by a factor of ~2 (SIO). By comparison, values of v^o are 
typically uncertain by ~10% - 20%. More accurate temperature discrimination is needed 
to search for any correlation of maximum plasma temperature with wind speed. High- 
resolution X-ray grating spectra can provide more detailed information on the distribution 
of X-ray plasma emission measure as a function of temperature. But improvements in X-ray 
telescope sensitivity will be needed to acquire grating spectra for a sample of WN stars that 
is sufficiently large to undertake correlation studies. 



5. Summary 

The main results of this study are the following: 

1. We report new X-ray detections of the WNL stars WR 16 (WN8h) and WR 78 (WN7h). 
These new detections, when combined with previous detections show that X-ray emis- 
sion is present in all WN subtypes from WN2 - WN9. This contrasts sharply with 
single carbon-rich WC stars which so far remain undetected in X-rays. 

2. The XMM detection of WR 16 is of low significance and the emission is too faint for 
spectral analysis. However, the X-ray spectrum of WR 78 reveals a high-temperature 
plasma component {kT hot ~ 2 keV), similar to that observed in other WN stars. Such 
hotter plasma is not predicted for X-ray emission arising solely in line-driven instability 
shocks. The origin of the hot plasma in WN stars is not yet known. Shock models 
predict that the maximum plasma temperature will scale with terminal wind speed 
according to kT^™^^^ oc v^. Such a dependence is not yet seen in the existing data 
but could be masked because temperatures determined from CCD X-ray spectra are 
average (not maximum) values and are in many cases uncertain by a factor of ~2. 

3. The X-ray luminosities of the WN8 stars WR 16 (a faint detection) and WR 40 (a 
non-detection) are an order of magnitude less than typical values observed for other 
detected WN stars. These WN8 stars also have the lowest terminal wind speeds of any 
WN stars observed so far, suggesting a possible link between terminal wind speed (or 
a parameter that depends on wind speed) and X-ray luminosity in single WN stars. 

4. The two new X-ray detections of WR 16 and WR 78 conform to and strengthen a 
previously noted trend of increasing Lx with wind luminosity Lwind = |Mv^ in WN 
stars. A regression fit based on existing data for 11 WN stars gives L^ oc L^^^^ °'*^ and 



-13- 



a correlation test yields a correlation probability Pcorr = 0.93. This trend needs to be 
confirmed in a larger sample of WN stars but tentatively suggests that wind kinetic 
energy may be an important factor in governing the X-ray luminosity levels of WN 

stars. 
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financial support from Bulgarian National Science Fund grant DO-02-85. 
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Table 1. XMM- Newton Observations (WNL Stars) 



Parameter 


WR 16 


WR 78 


Date 


2009 Dec. 28-29 


2010 Mar. 2-3 


Obsid 


0602020301 


0602020201 


Start Time (UTC) 


19:57:31 


14:08:02 


Stop Time (UTC) 


06:35:01 


00:45:08 


Optical Filter 


medium 


thick 


pn exposure (ks)^ 


30.0 


26.3 


MOS exposure (ks)^ 


33.0 


26.5 



''Usable exposure after removing high-background time 
intervals. MOS exposure is per MOS. The total exposure 
times before removing high-background intervals were 37.4 
ks (WR 16), 37.5 ks (WR 78). 
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Table 2. WNL Star Properties^ 



Name Sp. Type d Ay Voo log M log L^tar 

(kpc) (mag) (km/s) (Mo/yr) (Lq) 



WR 16 


WN8h 


3.60 


(2.37 - 


4.36) 


1.67 


670 


(630- 


740) 


-4.35 (- 


-4.55 - 


-4.20)*^ 


5.8 


(5.50 - 


6.15) 


WR78 


WN7h 


2.00 


(1.58- 


2.00) 


1.55 


1375 


(1365 ■ 


- 1385) 


-4.25 (- 


-4.42 - 


-4.10)<= 


5.9 


(5.63 - 


6.20) 


WR79a 


WN9ha 




1.99 




1.28 




955 




-4.60 (- 


-4.57 - 


-4.62)d 




5.78 





^ Spectral types are from van der Hucht (2001). Data for WR 16 and WR 78 show the values adopted in this study 
followed in parentheses by the range of values quoted in the literature (Crowther et al. 1995; Nugis & Lamers 2000; van 
der Hucht 2001; Hamann et al. 2006). Data for WR 79a are from Crowther & Bohannan (1997). 

'^ Clumping-corrected values are log M = —4.55 at d = 3.93 kpc from radio data (Nugis & Lamers 2000) and log M = 
—4.3 from empirical relations (Hamann et al. 2006). 

'^ Clumping-corrcctcd values are log M = —4.42 at d = 1.58 kpc from radio data (Nugis & Lamers 2000) and log M = 
—4.2 from empirical relations (Hamann ct al. 2006). 

•^Radio-derived and spectroscopic values are in good agreement, and clumping effects are thought to be small in WR 
79a (Crowther & Bohannan 1997). 
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Table 3. WNL Star X-ray Source Properties'' 



Name 



R.A. 
(J2000) 



Decl. 
(J2000) 



Counts'' 
(cts) 



Rate 

(cts/ks) 



Eso 
(kcV) 



log La;*^ Identification(oflfset) 

(ergs/s) (arcsec) 



WR 16 09 54 52.94 
WR 78 16 52 19.18 
WR 79a° 16 54 58.52 



57 43 39.2 


31±6 


1.03±0.20 


d 


d 


31.14 


41 51 15.8 


800±28 


30.5±1.10 


2.47 


0.12 (0.97) 


32.84 


41 09 02.9 


805±29 


23.0±0.83 


1.35 


0.63 (0.33) 


32.14 



GSC J095452. 90-574338.27 (0.96) 
GSC J165219. 25-415116.25 (0.93) 
GSC J165458. 50-410903.08 (0.30) 



^ Notes: XMM-Newton data are from EPIC pn. Parameters were determined using events in the 0.3 - 8 keV range. Tabulated quantities 
are: source name, J2000.0 X-ray position (R.A., Deck), pn net counts and net counts error (background-subtracted), pn count rate (Rate) 
obtained by dividing net counts by the usable exposure times given below, median photon energy E50, probability of constant count-rate 
(Peons) from the KS test followed in parenthesc by the value determined from x^ analysis of background-subtracted binned light curves 
using 1000 s bins, unabsorbed X-ray luminosity hx (0.3 - 8 keV) assuming the distances in Table 2, and HST GSC v2.3.2 optical counterpart 
identification. The offset (in arcsecs) between the X-ray and optical counterpart position is given in parentheses. The values of E50 and 
Pconst were evaluated using events within an extraction circle of radius 15" (68% encircled energy) in order to minimize background effects. 
Usable pn exposure times after removing high-background intervals were: 30.0 ks (WR 16), 26.26 ks (WR 78), and 35.03 ks (WR 79a). 

''Net counts for WR 78 and WR 79a are based on background-subtracted events within an extraction circle of radius 45" (90% encircled 
energy). A smaller extraction region of radius 15" (68% encircled energy) was used for the faint source WR 16 to minimize background 
contribution. 

"^Unabsorbed L^ values for WR 78 and WR 79a are based on fluxes measured in spectral fits after resetting the column density to zero 
(Njj = 0) for all model components. The value for WR 16 is based on a PIMMS simulation (Sec. 3.1). 

''insufficient counts for a reliable measurement. 

°Sce Skinner et al (2010) for additional information on WR 79a. 
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Table 4. X-ray Spectral Fits of WR 78 



Parameter 


Value 


Model 


NH,i*kTi + NH,2*kT2 


Abundances^ 


WN 


Nh,i (1022 cm-2) 


1.38 [1.20 - 1.60] 


kTi (keV) 


0.64 [0.48 - 0.88] 


normi (10"^) 


3.59 [2.40 - 6.30] 


Nh,2 (1022 cm-2) 


6.55 [4.10- 10.7] 


kTa (keV) 


2.27 [1.63 - 2.90] 


norm2 (lO^"^) 


14.5 [8.30 - 34.4] 


xVdof 


49.2/51 


xl 


0.96 


Fx (10"" ergs cm-2 g-i) 


2.15 (14.5) 


Fx,i (10-13 ergs cm-2 g-i-j 


0.33 (4.27) 


log Lx (ergs s"!) 


32.84 


log [Lx/Lstar] 


-6.64 


log [Lx/L„i„d] 


-4.68 



Note. — Based on simultaneous XSPEC (vers. 12.4.0) 
fit of the background-subtracted pn, MOSl, MOS2 spec- 
tra binned to a minimum of 20 counts per bin using the 
optically thin thermal plasma vapec model. All models 
included the XSPEC wabs photoelectric absorption com- 
ponent, based on Morrison & McCammon (1983) cross- 
sections and Anders & Ebihara (1982) relative abundances. 
The tabulated parameters are total equivalent neutral H 
absorption column density (Njj) including both ISM and 
local stellar (e.g. wind) contributions, the product of 
Boltzmann's constant time plasma temperature (kT), and 
XSPEC component normalization (norm) . Square brackets 
enclose 90% confidence intervals. The total X-ray flux (Fx) 
and fiux of the cool component (Fx.i) arc the absorbed 
values in the 0.3 - 8 keV range, followed in parentheses by 
unabsorbed values. The unabsorbed luminosity Lx (0.3 - 
8 keV) assumes d = 1.99 kpc (Table 2). hstar and hwind 
= (l/2)Mv^ are from Table 2. 

^Abundances were held fixed at the generic WN values 
given in Table 1 of Van der Hucht et al. (1986). The generic 
WN abundances reflect H depletion and N enrichment and 
are by number: He/H = 14.9, C/H = 1.90E-03, N/H = 
9.36E-02, O/H = 4.35E-03, Ne/H = 9.78E-03, Mg/H = 
3.26E-03, Si/H = 3.22E-03, P/H = 1.57E-05, S/H = 7.60E- 
04, Fe/H = 1.90E-03. All other elements were held fixed at 
solar abundances (Anders & Grevesse 1989). 
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Fig. 1.— Summed EPIC MOSl+2 contour image of WR 16 (0.3 - 8 keV; 33 ks of usable 
exposure per MOS; Gaussian-smoothed). The inner contour {2a) encloses a weak excess at 
the HST GSC stellar position (cross). 
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Fig. 2.— EPIC pn image of WR 78 (0.3 - 8 keV; 26.3 ks of usable exposure; 4."3 pixels) 
showing a clear detection at the HST GSC stellar position (cross). 
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Fig. 3. — Top and Middle: EPIC pn background-subtracted spectra of WR 78 and WR 79a, 
binned to a minimum of 5 counts per bin. Spectral fit parameters for WR 78 are given in 
Table 4. Best-fit parameters for WR 79a are Nh = 7.5 [6.4 - 8.3] x 10^^ cm'^ kTi = 0.55 
[0.46 - 0.61] keV, and kT2 = 2.66 [1.66 - 5.60] keV where brackets enclose 90% confidence 
ranges (SIO). Bottom: Unfolded spectrum of WR 78 binned to a minimum of 20 counts 
per bin showing the separate contributions of the cool (dotted line) and hot (dashed line) 
components (Table 4). The hot component accounts for most of the emission above 2 keV 
including the high-temperature Fe K line. Error bars omitted for clarity; log-log scale. 
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Fig. 4. — Unabsorbed X-ray luminosity (0.3 - 8 keV) versus stellar luminosity for WN 
stars. Assumed distances are given in Table 2 and SIO. WR star numbers are from van der 
Huclit 2001 and :li following the number denotes the presence of hydrogen in the spectrum. 
WNE and WNL classifications follow Hamann et al. 2006. Unabsorbed L^ was obtained by 
setting the column density to zero (Nh = 0) in best-fit spectral models and measuring the 
unabsorbed flux. Exceptions are the faint source WR 16 (PIMMS simulations; Sec. 3.1) and 
the undetected star WR 40. Unabsorbed Lj. values are from this work (WR 16, WR 78), 
Skinner et al. 2002 (WR 110), and Skinner et al. 2010 (and references therein). The L^ 
value for WR 1 is based on analysis of XMM archive data (obsid 0552220101) and assumes 
d = 1.82 kpc (van der Hucht 2001). The L^. upper limit for WR 40 is from Gosset et al. 
2005, adjusted upward by 0.14 dex to allow for distance uncertainties. Lgtar uncertainties 
reflect the range of values published in the literature (Crowther et al. 1995; Crowther & 
Bohannan 1997; Nugis & Lamers 2000; Hamann et al. 2006). The dashed line shows the 
canonical L^ = 10~'^Lf,oi relation for O and early B-type stars. The spectral types from van 
der Hucht (2001) are: WR 1 (WN4), WR 2 (WN2), WR 16 (WN8h), WR 18 (WN4), WR 
20b (WN6:h), WR 24 (WN6ha), WR 40 (WN8h), WR 78 (WN7h), WR 79a (WN9ha), WR 
110 (WN5), WR 134 (WN6), WR 136 (WN6h). 
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Fig. 5. — Unabsorbed X-ray luminosity (0.3 - 8 keV) versus wind luminosity Lu,ind = |Mv^ 
for WN stars. Mass-loss data are based on published values (see references below) and 
uncertainties reflect the range of pubhshed values. WR 20b is excluded for lack of reliable 
mass-loss data. The dashed line shows a linear regression fit log L^ = 0.89(± 0.47)log Ly^ind 
+ 26.41(±0.55) ergs s~^. Symbols and spectral types are the same as given in Figure 4. 
Mass-loss parameter references (in parentheses): WR 1: (4,7,10) WR 2: (1,2,4,8) WR 
16: (3,4,10); WR 18: (4); WR 24: (3,4,10); WR 40: (3,4,9,10) WR 78: (3,4,10); WR 79a: (5); 
WR 110: (4); WR 134: (1,4,10); WR 136: (1,2,4,6,7,10). References: (1) Abbott et al. 1986; 
(2) Howarth & Schmutz 1992; (3) Crowther et al. 1995; (4) Hamann et al. 2006; (5) 
Crowther & Bohannan 1997; (6) Ignace et al. 2001; (7) Ignace et al. 2003b; (8) van der 
Hucht 2001; (9) Gosset et al. 2005; (10) Nugis & Lamers 2000. 



